Weyl theory and explicit solutions of direct 
and inverse problems for a Dirac system with 
rectangular matrix potential 

B. Fritzsche, B. Kirstein, I.Ya. Roitberg, A.L. Sakhnovich 

Abstract 

A non-classical Weyl theory is developed for Dirac systems with 
rectangular matrix potentials. The notion of the Weyl function is 
introduced and the corresponding direct problem is treated. Further- 
more, explicit solutions of the direct and inverse problems are obtained 
for the case of rational Weyl matrix functions. 
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1 Introduction 

Consider self-adjoint Dirac-type (also called Dirac, ZS or AKNS) system, 
which is a classical matrix differential equation: 

-^y{x,z) = i{zj + jV{x))y{x,z) (x > 0), (1.1) 

(1.2) 



I mi 





, v = 


' 


V 







. ^* 


_ 



1 



where 1^^. is the x ruk identity matrix and v{x) is an mi x 7712 matrix 
function, which is called the potential of system. Dirac-type systems are very 
well-known in mathematics and applications (see, for instance, books [6, 7, 
30,32,47], recent pubhcations [3-5,9,10,17,18,48], and numerous references 
therein). The name ZS-AKNS is caused by the fact that system (1.1) is 
an auxiliary linear system for many important nonlinear integrable wave 
equations and as such it was studied, for instance, in [1,2, 13,21,42,52]. The 
Weyl and spectral theory of self-adjoint Dirac systems, where mi = 1712, was 
treated, for instance, in [5, 9, 24, 28, 32, 41, 47] (see also various references 
therein). The "non-classical" Weyl theory for the equally important case 
mi 7^ m2 and related questions are the subject of this paper. 

In Section 2 we treat the direct problem for the general-type Dirac sys- 
tem, that is, system (1.1), where the potential v is locally summable. A 
definition of the non-cxpansivc generalized Weyl function is given, its exis- 
tence and uniqueness arc proved, and some basic properties are studied. 

In Section 3 we consider Dirac systems with so called generalized pseudo- 
exponential potentials (see Definition 3.1). Direct and inverse problems for 
such systems are solved there explicitly. 

As usual, N stands for the set of natural numbers, R stands for the real 
axis, C stands for the complex plain, and C_|_ for the open upper semi-plane. 
If a G C, then a is its complex conjugate. The notation Im is used for 
image. An 1712 x mi matrix a is said to be non-expansive, if a*a < (or, 
equivalently, if aa* < 1^2)- 

We put mi -|- m2 =: m. The fundamental solution of system (1.1) is 
denoted by u{x, z), and this solution is normalized by the condition 

u{0,z)^lm. (1.3) 

2 Direct problem 

We consider Dirac system (1.1) on the semi-axis x G [0, 00) and assume 
that V is measurable and locally summable, that is, summable on all the 
finite intervals. In a way, which is similar, for instance, to the non-classical 
problem treated in [40] we shall use Mobius transformations and matrix balls 
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to solve the direct problem for Dirac system. 

Introduce a class of nonsingular m x mi matrix functions V{z) with 
property-j, which are an immediate analog of the classical pairs of parameter 
matrix functions. Namely, the matrix functions V{z) are meromorphic in 
C+ and satisfy (excluding, possibly, a discrete set of points) the following 
relations 

V{zYV{z) > 0, V{z)*3V{z) > {ze C+). (2.1) 

Definition 2.1 The setAf{x,z) of Mobius transformations is the set of val- 
ues at X, z of matrix functions 

^(x,z,V)^[0 Im,]u(x,z)-'V(z)(^[lm, o] u(x, zy^Viz))'', (2.2) 

where V{z) are nonsingular matrix functions with property-j. 

Proposition 2.2 Let Dirac system (1.1) on [0, oo) be given and assume 
that V is locally summable. Then the sets J\f{x, z) are well-defined. There is 
a unique matrix function (p{z) in C+ such that 

if{z)= f]Af{x,z). (2.3) 

x<oo 

This function is analytic and non-expansive. 
Proof. It is immediate from (1.1) that 

— (u{x, z)*ju{x, z)) = i{z — z)u{x, z)*u{x, z) < 0, z E C+. (2.4) 

According to (1.3) and (2.4) we have 

2l(x, z) = {Qlij{x, 2;)}-j=i := u{x, z)*ju{x, z) < j, z E C+, (2.5) 

where 21 is partitioned into four blocks so that Sljj is an mi x rrii matrix 
function {i — 1,2). Inequality (2.5) yields 

{u{x,z)*) ^ju{x,zy^ > j. (2.6) 
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Thus, we get 

det([V 0]u{x,z)-'V{z)') ^0, (2.7) 

and so Af is well-defined via (2.2). Indeed, if (2.7) does not hold, there is a 
vector / e C"*i such that 

[V 0] ju{x, z)-'V{z)f = [V 0] u{x, zr'V{z)f = 0, f^O. (2.8) 

By (2.1) and (2.6) the subspace Im {u{x, z)~^V{z)) is a maximal j-non- 
negative subspace. Clearly Im ( [/^^ O] ) is a maximal j-nonnegativc sub- 
space too. Therefore (2.8) implies u{x, z)~^V{z)f e Im ( [l^^ O]*)- 
then it follows from 

[I^, 0]uix,z)-'Viz)f = 

that / = 0, which contradicts the inequality in (2.8). 
Next, rewrite (2.2) in the equivalent form 

= u{x, z)-'V{z) ( [V 0] u{x, z)-'r{z)'^ ~\ (2.9) 

In view of (2.1), (2.9), and of the definition of 21 in (2.5), formula 

(p{z) eAf{x,z) (2.10) 

is equivalent to 

[V ^(^)*]2t(x,0) 

In a standard way, using formula (2.4) and the equivalence of (2.10) and 
(2.11), we get 

J\f{xi,z) gN'{x2,z) for xi>X2. (2-12) 
Moreover, (2.11) at x = means that 

X{0, z) = {(p{z) : ^{zy(p{z) < V}. (2.13) 



-'mi 

ifix.z.V)^ 



-'mi 



> 0. 



(2.11) 
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By Montel's theorem, formulas (2.12) and (2.13) imply the existence of an 
analytic and non-expansive matrix function (f{z) such that 



(2.14) 



x<oo 



Indeed, because of (2.12) and (2.13) we see that the set of functions ip{x, z, V) 
of the form (2.2) is uniformly bounded in C+. So, Montel's theorem is ap- 
plicable and there is an analytic matrix function, which we denote by <fioo{z) 
and which is a uniform limit of some sequence 



(Poo{z) = lim ip{xi, z, Vi) {i e N, Xi t, lim Xi = oo) 



(2.15) 



on all the bounded and closed subsets of C+. Since Xi t and equalities (2.9) 
and (2.12) hold, it follows that the matrix functions 



Vij{z) := u{xi,z) 



-'mi 



U > i) 



satisfy relations (2.1). Therefore, using (2.15) we derive that (2.1) holds for 



^i,oo(^) := u{xi,z) 



mi 



which, in its turn, implies that 



(2.16) 



Since (2.16) holds for all i e N, we see that (2.14) is true for ip^z) = (poo{z)- 
Now, let us show that jV is a matrix ball. It follows from (2.5) that 
2I22 < -Im2- Moreover, it follows from (2.4) and (2.5) that < 1(^-^)21 < 
i{z — z)j. Taking into account the inequalities above, we derive 

-Q122{X,Z) > {l + i{z- Z)x)lrr 



'm2 ■ 



(2.17) 



Note also that (2.5) implies Ql{x,z) ^ > j ior z e C+ (see [36]). Thus, we 
get 



(21-1)^^ = (2ln - 2li22l22^2l2i) ' > I„ 



(2.18) 
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As — 2I22 > the square root T = ( — 2I22) is well-defined and we rewrite 
(2.11) in the form 

- 2li22l2-2'2l2i - (^*T - 2li2T^i) (T^ - T-'^^i) > 0, 

where 2ti2 = ^21- Equivalently, we have 

(p = plUJpr - 21^2^2121, 00* U < /m2, (2-19) 

pi := = ( - 2122)"'^', Pr (2lii - 2li22l22^2l2i)'^'. (2.20) 

Here 00 is an m2 x mi matrix function. Since (2.10) is equivalent to (2.19), 
the sets A/'(x, z) (where the values of x and z are fixed) are matrix balls, 
indeed. According to (2.17), (2.18), and (2.20) the next formula holds: 

Pl{x,z)^0 (X^OO), Pr{x,z)<Imi- (2-21) 

Finally, relations (2.14), (2.19), and (2.21) imply (2.3). 



In view of Proposition 2.2 we define the Weyl function similar to the canonical 
systems case [47]. 

Definition 2.3 The Weyl-Titchmarsh (or simply Weyl) function of Dirac 
system (1.1) on [0, oo), where potential v is locally summable, is the function 
ip given by (2.3). 

By Proposition 2.2 the Weyl-Titchmarsh function always exists. Clearly, it 
is unique. 

Corollary 2.4 Let the conditions of Proposition 2.2 hold. Then the Weyl 
function is the unique function, which satisfies the inequality 



[Imi <p{zy]u{x,z)*u{x,z) 



_ip{z) 



dx < oo. 



(2.22) 



Proof. According to the equalities in (2.4) and (2.5) and to the inequality 
(2.11) we derive 



/ [Imi ^p{z)*]u{x,z)*u{x,z) 
Jo 



-'mi 



dx 



¥>{zr] (2.23) 



z — z 



X (21(0, z) - 2l(r, z)) 



-'mi 



< 



'—[Imi 'fi{z)*]^{0,z) 



z — z 



-'mi 
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Inequality (2.22) is immediate from (2.23). Moreover, as u*u > —21, the 
inequality (2.17) yields 



[0 Im,] 



U[X, Z) U[X, Z) 







'm2. 



dx > rlm2- 



(2.24) 



By (2.24), the function satisfying (2.22) is unique. 



Remark 2.5 In view of Corollary 2.4, inequality (2.22) can be used as an 
equivalent definition of the Weyl function. Definition of the form (2.22) is a 
more classical one and deals with solutions of (1.1) which belong to L^(0, oo). 
Compare with definitions of Weyl-Titchmarsh or M -functions for discrete 
and continuous systems in [10, 32, 33, 37, 38, ^1, 49, 50] (see also references 
therein). 

Our last proposition in this section is dedicated to a property of the Weyl 
function, the analog of which may be used as a definition of generalized Weyl 
functions in more complicated non-self-adjoint cases (see, e.g., [17,38,40]). 

Proposition 2.6 Let Dirac system (1.1) on [0, oo) be given, and assume 
that V is locally summable. Then, the following inequality 



sup 

x<i,zec+ 



^u{x, z) 



< oo 



(2.25) 



holds on any finite interval [0, I] for the Weyl function (p of this system. 

Proof. We fix some I. Now, choose x such that < a; < Z < oo. 
Because of (2.3) the Weyl function </? admits representations (2.2) (i.e., 
(p{z) — (f{x,z,V)). Hence, we can use (2.1) and (2.9) to get 



*(x, 2;)*j*(x, z) > 0, *(a:, z) := e"'^^u(a:, z) 



(2.26) 



On the other hand, equation (1.1) and definition of ^ in (2.26) imply that 
^ (e-2-^*(x, zYilm + mix, z)) = e-2-^*(x, zY 

X + ])]V - V3{Im + j)) - 2M{I^+ 3)Y{x,z) (2.27) 

'-2MI 



= 2e 



-■2xM 



^{x,zy 



'■nil 



IViX] 



-w{x)* 



'^(x,z), M := sup\\V(x) 

x<l 
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Using (2.26) and (2.27) we derive 
d 



4- (e-^^^^(x, zYil^ + j)^(x, z)] < 2e-2"^*(a;, zY 

•M/^^ ^(a;,^) < 0. (2.28) 
Finally, inequalities (2.26) and (2.28) lead us to 



X 



-w{xy 



^(x, zy^{x, z) < ^{x, zYilm + 3)^{x, z) < 2e2"^V, (2.29) 
and (2.25) follows. ■ 

3 Direct and inverse problem: explicit solu- 
tions 

Various versions of Backlund-Darboux transformations are actively used to 
construct explicit solutions of linear and intcgrablc nonlinear equations (see, 
e.g., [8,20,25,34,35,43,51] and numerous references therein). For the spectral 
and scattering results that follow from Backlund-Darboux transformations 
and related commutation and factorization methods see, for instance, pub- 
hcations [11, 12, 16, 19, 22, 24, 27, 29, 42]. Here we will give exphcit solutions 
of our direct and inverse problems using the GBDT version of the Backlund- 
Darboux transformation (see [14,16,24,39,42,43] and references therein). 
To obtain explicit solutions, we consider mi x m2 potentials v of the form 

v{x) = -2ii?i*e^^"*E(a;)-^e^^"i?2, (3.1) 

where some n G N is fixed and the n x n matrix function S is given by the 
formula 

E(x)=Eo+/ A{t)A{tydt (Eo>0), A(a;) - [e"^^"-!?! e^^°t?2] ■ (3.2) 

Here a, i^i, and are n x n, n x rrii, and n x m2 parameter matrices, and 
the following matrix identity holds: 

aSo - Soa* = i(i9i^?* - 192^^). (3.3) 

Clearly, E(a:) is invertible for x > and the potential v in (3.1) is well-defined. 
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Definition 3.1 The mi x 777,2 potentials v of the form (3.1), where relations 
(3.2) and (3.3) hold, are called generalized pseudo- exponential potentials. It 
is said that v is generated by the parameter matrices a, Eq, and '&2- 

According to [39, Theorem 3] (see also [16]), the fundamental solution u of 
system (1.1), where V is given by (1.2), v is a generalized pseudo-exponential 
potential, and u is normalized by (1.3), admits representation 

u{x, z) = Wa{x, z)e'''''^Wa{0, z)~^ . (3.4) 

Here we have 

Wa{x, z) := Im + ijA{xyE{x)-\zIn - a)-^A{x). (3.5) 

Note that the case mi = m2 (i.e., the case of the pseudo-exponential 
potentials) was treated in greater detail in [24] (see [24] and references therein 
for the term pseudo-exponential, itself, too). Formulas (3.2) and (3.3) yield 

aS(x) - J:{x)a* = iA{x)jA{x)*. (3.6) 

Identity (3.6), in its turn, implies that Wa{z) is a transfer matrix function 
in Lev Sakhnovich form [44-47]. However, Wa{x,z) possesses an additional 
variable x and the way, in which this matrix function depends on x, is essen- 
tial. 

From [16, formula (2.9)], where Wu and W21 are mi x mi and m2 x mi 
blocks of Wa{0, z): 

w^{0,z)^:{W,j{z)}l^„ (3.7) 

we see that 

W2i{z)Wn{z)-^ ^ -i^;j:o\zin - ey^-du e -.^ a - id^-^in-K (3.8) 

Wc note that [16, formulas (2.6) and (2.7)] imply that Wii{z) is always well- 
defined and invertible for z ^ a{a) U ^{O), where a denotes spectrum. 

Relations (3.4), (3.7), and (3.8) are basic to solve the direct problem for 
Dirac systems with the generalized pseudo-exponential potentials (3.1). 
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Theorem 3.2 Let Dirac system (1.1) on [0, oo) be given and assume that 

V is a generalized pseudo-exponential potential, which is generated by the 
matrices a, Eq, and '§2- Then the Weyl function ip of system (1.1) has 
the form: 

ip{z) = -i'd;E^\zin - ey^^i, e ^ a - i-dirji^ K (3.9) 

Proof. We compare (3.8) and (3.9) to see that 

^{z) = W2i{z)Wi^{z)-\ (3.10) 
Because of (3.4), (3.7), and (3.10) we have 

W^^{z)-\ (3.11) 

To consider the matrix function A*E~^, which appears in the definition (3.5) 
of Wai we derive from (3.2) that 

S(x)-^A(x)A(x)*S(x)-^ = -^T.{x)-\ (3.12) 
It is immediate also from (3.2) that E(x) > 0. Therefore, using (3.12) we get 

poo 

/ E(t)-^A(t)A(t)*E(i)-Mt < Eo (3.13) 
Jo 

Furthermore, the last equality in (3.2) implies that 

sup ||e^^"A(x)|| < M, (£>0). (3.14) 

9z>||a||+£ 

It follows from (3.5), (3.13), and (3.14) that the entries of the right-hand side 
of (3.11) are well-defined and uniformly bounded in the iv^(0, oo) norm with 
respect to x for all z such that '^z > max (||q;||, ||^||) -|- e and £ > 0. Hence, 
taking into account (3.11) we see that (2.22) holds for z from the mentioned 
above domain. So, according to the uniqueness statement in Corollary 2.4, 
ip{z) of the form (3.9) coincides with the Weyl function in that domain. Since 



I mi 

^{z)_ 



e'''^Wa{x,z) 



'■nil 
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the Weyl function is analytic in C+, the matrix function coincides with it 
in C+ (i.e., is the Weyl function, indeed). ■ 

For the case that v is a generalized pseudo-exponential potential, where Eq > 
0, our Weyl function coincides with the reflection coefficient from [16] (see [16, 
Theorem 3.3]). Hence, the solution of our inverse problem can be considered 
as a particular case of the solution of the inverse problem from [16, Theorem 
4.1], where the singular case Eq was treated too. 

Before we formulate the procedure to solve inverse problem, some results 
on rational matrix functions and notions from system and control theories 
are required (see, e.g., [26,31]). Let (f{z) be a strictly proper rational matrix 
function, that is, such a rational matrix function that 

lim ^p{z) = 0. (3.15) 

Then p admits representations (also called realizations): 

^{z) = Cj,{zIm - AnY^Bj,, (3.16) 

where Cjv, An^ and B]^ are 1112 x N, N x N, and x mi, respectively, 
matrices. Here N E N, and mi (m2) denotes the number of columns (rows) 
of (f. 

Definition 3.3 The minimal possible value of N in realizations (3.16) is 
called the McMillan degree of ip, and we denote this value by n. Realizations 
(3.16), where N = n, are called the minimal realizations. 

Prom [31, Theorems 21.1.3, 21.2.1] we easily see that for a minimal realization 

ip{z) = C{zIn-A)-^B (3.17) 

of a matrix which is non-expansive on IR and has no poles in C+, there is 
a positive solution X > of the Riccati equation 

XC*CX + i{XA* - AX) + BB* = Q. (3.18) 

Furthermore, all the hermitian solutions of (3.18) are positive. 
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Theorem 3.4 Let (f{z) be a strictly proper rational matrix function, which 

is non- expansive on M. and has no poles in C+. Assume that (3.17) is its 
minimal realization and that X > is a solution of (3.18). 

Then ip{z) is the Weyl function of Dirac system, the potential of which 
is given by (3.1) and (3.2), where 



A + iBB*X-^, Sq = X, = B, = -i^C*. (3.19) 



This solution of the inverse problem is unique in the class of Dirac systems 
with locally hounded potentials. 

Proof. From (3.19) we see that 

aEo - Soa* ^ AX - XA* + 2iBB\ - i?2^;) = iBB* - \XC*CX, 

and so (3.3) is equivalent to (3.18). Since (3.3) holds, we apply Theorem 
3.2. Theorem 3.2 implies that the Weyl function of the Dirac system, where 
V is given by (3.1), has the form (3.9). Next, we substitute (3.19) into 
(3.9), to derive that the right-hand sides of (3.17) and the first equality 
in (3.9) coincide. In other words, the Weyl function of our system admits 
representation (3.17). 

Finally, the uniqueness of the solution of the inverse problem follows from 
Theorem [15, Theorem 4.1]. ■ 

We note that the corresponding uniqueness result in [16] was proved only for 
the class of systems with the generalized pseudo-exponential potentials. 

Because of the second equality in (3.8) and identity (3.3), the matrix 6 
satisfies another identity: 6'Eo — Eq^* = —\{§i'd\ -\- •&2'^2)^ ^^^^ 



If / 7^ is an eigenvector of 6 (i.e.. Of — A/), identity (3.20) implies that 



En - rE, 



'0 ^ 



iEo-i(^?ii?t + ^?2^?;)E, 



'0 ^ 



(3.20) 



(A-A)rEoV 



irEo-^(i?i^?t+i?2^?;)Eo-V- 



(3.21) 



Since Eq > 0, we derive from (3.21) that 



a{e) C C_ U M. 



(3.22) 
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Real eigenvalues of 9 play a special role in the spectral theory of an 
operator, which corresponds to the Dirac system with a generalized pseudo- 
exponential potential (see, e.g., [23] for the case of square potentials). In our 
case the operator "H corresponding to the Dirac system is defined in a way, 
which is similar to the definition from [23], but the initial condition is quite 
different. Namely, we determine T-L by the differential expression 

n^^y = -ij^y - Vy, (3.23) 

and by its domain ViJ-t), which consists of all locally absolutely continuous 
C"- valued functions y in 1/^(0, oo), such that 

ndeyeLl{0,oo), y{0) = 0. (3.24) 

Proposition 3.5 Let the conditions of Theorem 3.2 hold, let 9 be given by 
the second relation in (3.8), and let A be an eigenvalue of 9: 

Of^Xf, f^O, XeR. (3.25) 

Then, the matrix function 

g{x):^jA{xyi:{x)-'f (3.26) 
is a bounded state ofH and T-Lg — Xg. 

Proof. First, we show that formulas (3.21) and (3.25) yield 

i?tSo-V = 0, r.E^'f^O, af^Xf. (3.27) 

Indeed, the first two equalities in (3.27) easily follow from (3.21) for the case 
that A = A. The equality af = Xf is immediate from 9f — A/, definition of 
9 in (3.8), and equality i^^Sq V = 0. 
Next, we show that 

(jA*E-^)' = i/A*E-^a + (A*E-^A - jA*E-^Aj)jA*E-^ (3.28) 

Formula (3.28) follows from a general GBDT formula [43, (3.14)] and also 
from its Dirac system subcase [43, (2.13)], but it will be convenient to prove 
(3.28) directly. We note that formula (3.2) imphes 

A' = -iaAj, E' = AA*, (3.29) 
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and formula (3.6) can be rewritten as 

a*T.-^ = E-^a - iE'^AjA^E-^ (3.30) 

Since = using (3.29) and (3.30) we obtain (3.28). 

Now, partitioning A into two blocks and using (3.1) and (3.2), we see 
that 

v(x) = -2iAi(x)*E(a;)"^A2(x), A =: [Ai Aa] . (3.31) 

In view of (1.2) and (3.31) we have 

A*E-^A - jA*E-^Aj = ijV. (3.32) 

Applying both sides of (3.28) to / and taking into account the last equahty 
in (3.27) and relation (3.32), we derive 

(jA(a;)*E(a;)-V)' = iXfA{xrE{x)-'f + ijV(a;)jA(x)*E(x)- V- (3.33) 

Because of (3.23) and (3.26), we can rewrite (3.33) as 

Hd.g = Xg, (3.34) 

and it remains to show that g e T>{H), that is, that g e L'^{0, oo) and (3.24) 
holds for y — g. Prom (3.13) and (3.34) we see that g, Hdeg £ L'^i^i 
Finally, the initial condition 

g{0)= [^i -t?2]*So-V = (3.35) 
is immediate from (3.2), (3.26), and (3.27). ■ 
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